First observation of 13 Li ground state 
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The ground state of neutron-rich unbound 13 Li was observed for the first time in a one-proton re- 
moval reaction from 14 Be at a beam energy of 53.6 MeV/u. The 13 Li ground state was reconstructed 
from 11 Li and two neutrons giving a resonance energy of 120lgQ keV. All events involving single 
and double neutron interactions in the Modular Neutron Array (MoNA) were analyzed, simulated, 
and fitted self-consistently. The three-body ( 11 Li+n+n) correlations within Jacobi coordinates 
showed strong dineutron characteristics. The decay energy spectrum of the intermediate 12 Li sys- 
tem ( 11 Li+n) was described with an s-wave scattering length of greater than —4 fm, which is a 
smaller absolute value than reported in a previous measurement. 

PACS numbers: 21.10.Dr, 25.60.-t, 29.30.Hs 



The increasing availability of rare isotope beams has 
made it possible to extend nuclear structure measure- 
ments to nuclei far away from stability. In light neutron- 
rich nuclei these studies have been extended beyond the 
neutron dripline [l|. The existence of two- neutron halos, 
first observed in 11 Li @, 0], sparked the interest in un- 
bound two-neutron systems where measurements of the 
decay products could yield information about neutron- 
neutron correlations inside the nucleus. Recent examples 
include 10 He @,[|, 13 Li 0,0, 16 Be 0, and 26 0]. In 
16 Be the first evidence for ground state dineutron decay 
was observed because sequential decay via 15 Be was en- 
ergetically not possible 0. 

In this article, we report on a new measurement of the 
13 Li decay products. Very few calculations exist that pre- 
dict the energy levels of 13 Li. Over 25 years ago, a shell 
model calculation predicted the ground state of 13 Li to 
have a spin and parity of § located about 3 MeV above 
the n Li+2?i threshold @. In 2008 Aksyutina et al. @ re- 
ported a resonance at a decay energy of 1.47(31) MeV in 
a proton removal reaction from 14 Be at 304 MeV/u with a 
liquid hydrogen target. In the same reaction the unbound 
intermediate nucleus 12 Li was observed to decay via the 
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emission of an s-wave neutron with a scattering length 
of —13.7(1.6) fm. We repeated this proton removal reac- 
tion at a 14 Be beam energy of 53.6 MeV/u on a beryllium 
target. In addition to measuring a new low-lying reso- 
nance in 13 Li at 120t§o keV, which is presumably the 
ground state, we extract a limit on the scattering length 
of a s > —4 fm for 12 Li. 

The experiment was performed at the Coupled Cy- 
clotron Facility of the National Superconducting Cy- 
clotron Laboratory (NSCL) at Michigan State Univer- 
sity. A 120 MeV/u 18 primary beam bombarded 
a 4113 mg/cm 2 Be production target to produce the 
53.6 MeV/u 14 Bc secondary beam. The A1900 fragment 
separator, with a 1050 mg/cm 2 Al achromatic wedge de- 
grader, was used to separate the 14 Be secondary beam 
from other reaction products and the primary beam. 
The rate of the 14 Be secondary beam was about 500 pps 
at the reaction target and had a momentum spread of 
2.5%. Contaminant nuclei in the secondary beam were 
excluded from analysis based on event-by-event measure- 
ment of the time-of-flight from plastic scintillator tim- 
ing detectors located at the A1900 focal plane and just 
before the reaction target. The 14 Bc beam impinged 
on a 477 mg/cm 2 Be reaction target and produced 13 Li 
through one-proton removal. The two-neutron unbound 
13 Li decayed immediately into u Li and two neutrons. 

The 11 Li fragments were deflected by the Sweeper su- 
perconducting dipole magnet fioj that was set to a mag- 
netic rigidity of 3.735 Tm. The detection setup was 
identical to the experiment by Hall et al. 
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measured excited states in 12 Li following the two-proton 
removal from 14 B. The deflected nuclei passed through 
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two x-y position sensitive cathode-readout drift chambers 
(CRDCs), separated by 1.816 m, a 5 mm thick plastic 
scintillator, and were stopped in a 150 mm thick plastic 
scintillator. An inverse ion-optical matrix for ray tracing 
was created with the program COSY infinity [12| using 
the measured magnetic field map of the Sweeper to pro- 
vide tracking of the n Li trajectories back to the target 
position. 

The neutrons produced from the 13 Li breakup were 
detected using the Modular Neutron Array (MoNA) fl3l . 
[lij . The individual plastic scintillator bars were arranged 
in nine walls, each 16 bars high, resulting in an active area 
2 m wide by 1.6 m tall. The distance from the target to 
the center of the first layer was 8.44 m. 

The three-body decay energy of 13 Li can be calculated 
with £ D ccay = Mi3 Li -Mii L i-2M n! where Mi 3Li (Mn Li ) 
is the mass of 13 Li ( 11 Li) and M n is the neutron mass. 
The invariant mass, Mi3Li, was calculated from the ex- 
perimentally measured four-momenta of the 11 Li and two 
neutrons. It should be noted that all experimental de- 
cay energy spectra are presented without unfolding of the 
resolution or acceptance. 

The reconstruction of the 13 Li decay energy relies 
on the correct identification of two neutrons interacting 
in MoNA. It is critical to distinguish real two-neutron 
events from events where a single neutron re-scatters and 
interacts twice in MoNA. It has been shown that causal- 
ity cuts applied to the distance and the velocity between 
two successive interactions can effectively enhance the se- 
lection of real two- neutron events @, H, EH, EH • Fig- Q] 
shows the reconstructed three-body decay energy spec- 
trum of 13 Li ( 11 Li+2n) without (a) and with (b) the 
causality cuts applied. The cuts require that two succes- 
sive interactions be separated by more than 50 cm and 
that the time between the first and second interaction be 
shorter than the time it would take a neutron at beam ve- 
locity to travel between the respective interaction points 
(see Ref. [l7| and references therein for additional infor- 
mation on In cuts). The data exhibit a resonance-like 
structure below 500 kcV, which was not observed in the 
previous measurement of 13 Li [f|. However, the earlier 
experiment had zero efficiency for the detection of two- 
neutron events below 200 keV [lH , which was also shown 
to have significant effects on the reconstructed shape of 
the decay energy spectrum of the excited 11 Li fl5| . 

In order to analyze the observed peak for the pres- 
ence of a resonance, detailed Monte Carlo simulations 
were performed 
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The simulations included the in- 
coming beam characteristics, the reaction mechanism to 
populate specific states, and their subsequent decay as 
well as the detector resolutions and efficiencies. The one- 
proton and pn-removal reactions were simulated using 
the Glauber reaction model. As mentioned earlier the 
detailed response of MoNA to the interactions of multi- 
ple neutrons is especially important and was simulated 
with Geant4 [l|| [2(| using the custom neutron interac- 
tion model menate_r [2l|. In contrast to the standard 
intranuclear cascade models in Geant4, which uses only 



tn 
+■» 

E 
3 

o 
o 



c 

3 
O 

o 



250 
200 
150 
100 
50 


60 
50 
40 
30 
20 
10 

0. 



i i "Li + 2n 


(a) 






f \ u 


• Exp. 


Sum 


,.- — \— 


,3 Li 

...... «m 


F ^ '""m, 

... T .... T .... jT.TTTS-T 




11 i - n 1-0 

: Li + 2n 
La Csusslitv Cuts 

^ A. 
h X 


ji (b) 


X 


1 2 3 

E decay (MeV) 











0.5 



1.5 



2.5 



3.5 



"Decay 



(MeV) 



FIG. 1. (color online) Measured three-body decay energy 
spectrum of Li+2n without (a) and with (b) causality cuts 
as described in the text. The experimental data is fit with two 
components: (1) the dineutron decay of 13 Li— > 11 Li+2n (red 
long-dashed line) and (2) the 12 Li— > 11 Li+n decay (blue dot- 
dashed line) The input 13 Li resonance distribution is shown 
in the insert of panel (b). 



total inelastic reaction cross sections for neutrons above 
20 MeV, menate_R uses cross sections for the differ- 
ent reaction channels, including (n,np), (n,p), (n, rvy), 
and (n, a). It has recently been shown that the use of 
MENATE_r is important for correctly simulating the re- 
sponse of plastic scintillator detectors 17[ • Distortions of 
the experimental decay energy spectra are present due to 
the efficiency and acceptance of the detector setup [see 
the inserts of Figs. HJb) and [Ha)]. The accuracy of the 
Monte Carlo simulation was estimated through varying 
input parameters, such as the momentum distribution 
of the residual 11 Li fragment, and examining how this 
changed the efficiency curves. The results showed that 
large variations of the input parameters resulted in mini- 
mal changes (less than 5%) to the shape of the efficiency 
curves. 



In addition to the three-body decay energy ( 13 Li), the 
simulations had to consistently describe other observ- 
ables, such as the decay of 12 Li, which could be directly 
populated from pn-stripping, and the overall multiplic- 
ity distribution of interactions in MoNA. Limiting the 
data to multiplicity = 1 events should enhance the con- 
tributions of the direct 12 Li population. Fig.[2]shows the 
overall two-body ( 11 Li+n,) decay energy spectrum (a), 
the multiplicity = 1 gated n Li+n decay energy spec- 
trum (b), as well as the multiplicity distribution (c) for 
events in coincidence with 11 Li. The two-body decay en- 
ergy was calculated as i^Decay = ^ 12 Li — — M n , 
where the invariant 12 Li mass was calculated from the 
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FIG. 2. (color online) Measured decay energy spectrum of Li+n including all multiplicities (a) and with a multiplicity = 1 
gate (b). Panel (c) displays the overall multiplicity distribution for all events in coincidence with a xl Li fragment. The insert 
in panel (a) shows the efficiency for detecting one neutron from a In decay (solid line) and two neutrons from a 2n decay with 
(dashed line) and without (long-dashed line) the causality cuts. The 2n decay was simulated with each neutron being emitted 
with 50% of the total decay energy. The \ 2 as a function of the 12 Li s-wave scattering length is shown in the insert of panel (b). 



detected 11 Li fragment and neutron. 

The Monte Carlo simulation included two components: 
(1) the decay of 13 Li from the one-proton removal reac- 
tion and (2) the decay of 12 Li from pn-stripping. The 
black solid lines in Figs. [1] and [2] represent the sum of the 
13 Li and 12 Li contributions to the corresponding spec- 
tra. As shown in Fig. (Ua), a component of the decay 
of 12 Li— > n Li+n (blue dot-dashed line) is present in the 
three-body decay energy spectrum from events in which 
the single neutron interacted twice in MoNA and thus 
added false two neutron events to the 13 Li spectrum. Ap- 
plying the causality cuts to Fig. QJa) should greatly re- 
duce the presence of false two neutron events in the three- 
body decay spectrum. The effectiveness of the causality 
cuts can be seen from the near absence of the 12 Li com- 
ponent in Fig. [ljb). Similarly, in the multiplicity = 1 
two-body decay spectrum, Fig.[2jb), a component of the 
13 Li— > n Li+2n decay (red dashed line) is present from 
events in which only one neutron is detected. 

The correlations of the three-body lx Li+n+n system 
can provide additional insight into the structure and de- 
cay of 13 Li. A complete description of the three-body 
correlations can be obtained from the relative energy 
(E x /Et) and the angle (9k) calculated within the T and 
Y Jacobi systems. The relative energy is defined as the 
energy of the two-body system (frag+n or n+n), E x , 
relative to the total three-body energy, Ey- A detailed 
description and illustration of the Jacobi coordinate sys- 
tems can be found in Refs. [22TE5I ] . In Fig. [3] the rel- 
ative energy and angle for both the T and Y systems 
are shown from the experimental data with the causal- 
ity cuts applied and compared to a three-body phase 
space and dineutron simulated emission. The three- 
body phase space distribution was calculated according 
to Refs. [1^, [27j with a Breit Wigner lineshape used to 
fit the three-body decay energy. Further details of the 
dineutron decay calculation are provided later in the text. 

The relative energy in the T system and angle in the 



Y system are particularly sensitive to the decay mecha- 
nisms, three-body phase space versus dineutron. In the T 
system E x /Et represents the energy of the n-n system 
relative to the total three-body energy. In the Y sys- 
tem 9k is the angle between the frag+n center-of-mass 
and the additional neutron. The difference between the 
three-body phase space and dineutron decay simulations 
can be clearly seen in the top-left and bottom-right pan- 
els of Fig. [3] While the phase space emission (red dashed 
line) does not fit the experimental data, the dineutron 
decay (green solid line) agrees well with the data show- 
ing a strongly correlated n — n emission. This presents 
the second case (the first being 16 Be Q) in which a 
ground state dineutron-like emission has been observed. 
A strong dineutron component (~ 50%) in the decay of 
13 Li for -Ebocay > 3 McV was also reported by Johansson 
et al. [|. 

Since the three-body correlations indicate that the 
dominate decay channel of 13 Li is the emission of a dineu- 
tron (or two strongly correlated neutrons), the two-body 
and three-body decay energy spectra need to be analyzed 
accordingly. Thus, the decay of 13 Li was simulated as the 
emission of a dineutron, while the 12 Li decay was simu- 
lated through the emission of a s-wave neutron with an 
associated scattering length. 

The dineutron decay was simulated as a two-step pro- 
cess following the formulism of Ref. [28| where the dineu- 
tron was emitted with energy E y (in the T, "cluster" sys- 
tem) from the 13 Li followed by the decay of the dineutron 
with energy E x ( 13 Li -> n Li + 2 n -> u Li + n + n). The 
total amplitude for the two-step dineutron decay process 
is constructed from the one-body decay amplitudes as a 
second order process [29| : 



A(E y ,E x ) 



A 1 (E y )A 2 (E x ) 
E x - (E v - ^T V {E X )Y 



(1) 



here Ai(E y ) and A 2 (E X ) represent the amplitudes of the 
dineutron emission with kinetic energy E y and the sub- 
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FIG. 3. (color online) The relative energy (E x /Et) and angle 
(cos(0fc)) as denned in the T and Y systems are shown for the 
experimental data (black filled circles) with causality cuts, 
three-body phase space decay simulation (red-dashed line), 
and dineutron decay simulation (green solid line). 



sequent dineutron breakup with intrinsic energy E x , re- 
spectively. Ey and Ty(E x ) represent the energy and 
width parameters of the propagator describing the dineu- 
tron virtual state. These parameters along with A2(E X ) 
are taken from the low-energy n — n scattering theory 
in free space, where the scattering length a s = —18.7 
fm (30j . assuming that the dineutron is a 1 sq n — n state. 

The Fermi Golden Rule gives the partial decay width 
distribution for the sequential process as 



^ = 2ttS{E t - E y - E x ) \A(E y ,E x )\ 2 



dEydE.. 



(2) 



where Ex, again, represents the total energy of the three- 
body system. The probability distribution then follows 
the usual Brcit-Wigncr form 



dP(E yi E x ) 
dEy dE x 



i 



dT{E T ) 

{E T - E r y + Tl{E T )/A dE y dE x 



(3) 



where E r and Y r (E) are the resonance energy and total 
width of the initial state in 13 Li. The delta-function in 
Eq. © enforces that E T = E y + E x in Eq. ©. From 
the presented formalism, the 13 Li decay was simulated as 
the emission of a dineutron with kinetic energy E y which 
proceeded to breakup with energy E x . 

The decay of 12 Li, which was populated through a 
pn— removal, was simulated usin g a s -wave lineshape cal- 
culated from the model of Ref. [3l|. The potential pa- 
rameters used in Ref. [3l| to describe the 14 Be and 14 B 
projectiles are used in the present work. This calculation 
was also used to describe the s-wave lineshape of the 12 Li 
decay populated from the 14 B(-2j?) reaction of Ref. (Tlj . 

The four decay energy spectra of Figs. Q] and [5] were 
fit simultaneously varying the 13 Li resonance parameters 
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FIG. 4. (color online) Measured decay energy spectrum of 
12 Li populated in the two-proton removal reaction from 14 B. 
Three components are shown to fit the data: two Breit- 
Wigner resonances (BW) and a —3.8 fm scattering length 
s-wave, corresponding to the extracted a 3 > —4 fm limit 
from the current work. Experimental data was taken from 
Ref. El. 



[E r and T r (E) in Eq. and the 12 Li s-wave scatter- 
ing length. The simultaneous fitting of all the three- and 
two-body decay spectra greatly increased the sensitiv- 
ity and robustness of the fit in comparison to individ- 
ual fits of each spectrum. The best overall fit (\ 2 min- 
imization) was achieved with a resonance state in 13 Li 
at E-l = 120±|g keV with T r = 125±% keV. A limit 
on the scattering length of a s > —4 fm for the 12 Li s- 
wave was determined from the % 2 fit which is shown in 
the insert of Fig. [2jb). A reduced chi-square (x 2 /^) of 
1.18 was obtained from the fit with 181 degrees of free- 
dom. The mass excess of 13 Li can be determined from 



the 120lgQ keV resonance and the mass excess of 11 Li 
(40.72828(64) MeV H) to be 56.99(8) MeV. 

The extracted scattering length limit of a s > —4 fm 
for 12 Li from the present data seems to be inconsistent 
with the previous result of —13.7(1.6) fm by Aksyutina 
et al. Q]. While Aksyutina et al. used the calculation 
presented in Ref. @ to calculate the s-wave lineshape in 
comparison to Ref. [3llj used in the present work, it was 
verified that the difference in the calculated lincshapes 
could not account for the difference in the extracted scat- 
tering lengths. However, one possible explanation for the 
difference would be that the events from the unidentified 
low-energy peak in 13 Li, due to the zero two-neutron effi- 
ciency below 200 keV, were included in the 12 Li spectrum 
of Ref. [6j . This would enhance the yield at low energy 
thus shifting the extracted scattering length of the fit to 
larger negative values. 

A previous experiment populating 12 Li in a two- 
proton removal reaction from 14 B included the -13.7 fm 
scattering length s-wave from Aksyutina et al. [6j along 
with two narrow resonances in the fit of the 12 Li decay 
energy spectrum. In order to have a consistent analysis, 
the 12 Li decay energy spectrum from Ref. [ll| has to be 
re-fit using an s-wave scattering length greater than —4 
fm, as extracted in the current work. Fig.5]shows the re- 
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FIG. 5. Level scheme of Li and Li determined from the 
current work. The 1470 keV resonance in 13 Li is from Ref. [(|. 
The broad s-state with a s > —4 fm in 12 Li is not shown. 



suits of the fit (x 2 jv — 1.08) where the s-wave was fixed 
at a scattering length of —3.8 fm, while the resonance pa- 
rameters of the other two previously observed resonances 
were free fit parameters. The values of the lower (pink 
long-dashed lines) and higher (green short-dashed line) 
resonances had to be adjusted slightly from 250(20) keV 
to 210(30) keV and from 555(20) keV to 525(25), respec- 
tively. 

One-neutron removal experiments and shell model cal- 
culations have indicated that the 14 B ground state has 
a dominant (64% — 89%) Si/ 2 configuration [3l]. I33T - I35] . 
Thus, a strong population of the 12 Li s-wave from 14 B 
would be expected along with a smaller population of 
d-w&ve resonances. While the fit shown in Fig. @] has 
a 50% s-wave component, the statistical uncertainty of 
the data allows for fits of similar quality (albeit having 
slightly larger \ 2 ) to be obtained with s-wave compo- 
nents ranging from 45% to 90%, which is in agreement 
with s and d components of the 14 B ground state. In 
populating 12 Li from 14 Bc(-pn), inclusion of the two nar- 
row resonances in the fit of Figs. Q] and [2] did not have a 
significant effect and slightly increased the x 2 jv to 1.2. 
Thus, the narrow resonances were not included in the 
presented fit of Figs. [T] and [2] The results suggest that 
14 Be must have a dominant S1/2 configuration, which is in 
agreement with the Coulomb dissociation measurements 
of Labiche et al. (3(| • 

The new level schemes of 13 Li and 12 Li, based on the 
present data and the data of Refs. @ and [ll|, are shown 
in Fig. [5] The present experiment has very low detection 
efficiency for two neutron events above 1 MeV and is thus 
insensitive to the previously observed level at 1.47 MeV 



in 13 Li [6J]. The 12 Li scattering length is not shown in 
Fig. [5] since it corresponds to a very broad distribution 
in excitation energy [lj. Thus the state at 210 keV has to 
be considered the ground state of 12 Li. The mass excess 
of 12 Li is then 49.01(3) MeV using the most recent mass 
measurements of 11 Li [32l |. 

The 210 keV ground state of 12 Li produces a scenario 
in which 13 Li is bound with respect to In emission and 
unbound with respect to 2n emission. The decay of the 
13 Li ground state is, thus, a genuine three-body process 
and, while a purely sequential decay through the n — n 
system is unlikely, the three-body correlations indicate 
the presence of the n — n virtual state (a» = -18.7 fm). 
The recent work by I. A. Egorova et al. [33] has demon- 
strated the complexity present in the three-body dynam- 
ics of the two-proton decay of 6 Be through the use of 
a three-cluster model and similar calculations could pro- 
vide new insights into the two-neutron emission from 13 Li 
and 16 Bc 0. 

In summary, the population and decay of 13 Li and 
12 Li were measured in p— and pn— removal reactions 
from 14 Bc. A low-lying resonance of 120^gQ keV was 
observed in 13 Li for the first time, corresponding to a 
mass excess of 56.99(8) MeV. Within a consistent de- 
scription of all measured decay spectra the decay of 
12 Li could be described by an s-wave with a scattering 
length of greater than —4 fm, which differs significantly 
from the —13.7(1.6) fm scattering length extracted in 
Ref. 0. The a s > —4 fm scattering length limit, from 
this work, was used to re-fit the decay energy spectrum of 
12 Li pop ulated with a two-proton removal reaction from 
14 B [11(. The observed 210 keV resonance in 12 Li repre- 
sents the ground state and corresponds to a mass excess 
of 49.01(3) MeV. This 12 Li ground state produces condi- 
tions favorable for true two-nucleon emission from 13 Li to 
n Li. The angular and energy correlations in the three- 
body system of 13 Li showed strong dineutron character 
for the ground state decay. 
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